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ABSTRACT 

Optical  lattices  are  employed  to  simulate  important  condensed  matter  materials,  with  ultracold  atoms  playing  the  role  of  the  electrons  in  a 
real  material.  We  have  obtained  the  phase  diagram  of  a  one-dimensional  (ID)  Fermi  gas  with  spin-imbalance  in  this  way.  This  system 
features  a  partially  polarized  phase  which  is  the  ID  analog  of  the  long-sought  FFLO  state.  Current  research  is  focused  on  direct  detection  of 
the  FFLO  state  and  explore  the  1D-3D  dimensional  crossover.  We  have  recently  detected  antiferromagnetic  order  in  the  3D  Fermi-FIubbard 
model  using  spin-dependent  Bragg  scattering  of  light.  This  was  a  major  goal  of  our  project,  enabled  by  the  development  of  a  novel  scheme 
to  evaporatively  cool  atoms  in  an  optical  lattice.  We  have  also  focused  on  developing  alternative  cooling  schemes.  We  have  also 
investigated  the  connection  between  quantum  integrability  and  thermalization  in  a  ID  Bose  gas,  and  considered  equilibration  and 
thermalization  phenomena  in  optical  lattices.  Other  theoretical  issues,  such  as  the  effects  of  spin-orbit  coupling,  the  disordered 
Bose-FIubbard  model,  and  many-body  localization,  have  also  been  explored. 
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Abstract 

Optical  lattices  are  employed  to  simulate  important  condensed  matter  materials,  with 
ultracold  atoms  playing  the  role  of  the  electrons  in  a  real  material.  We  have  obtained  the  phase 
diagram  of  a  one-dimensional  (ID)  Fermi  gas  with  spin-imbalance  in  this  way.  This  system 
features  a  partially  polarized  phase  which  is  the  ID  analog  of  the  long-sought  FFLO  state. 
Current  research  is  focused  on  direct  detection  of  the  FFLO  state  and  explore  the  1D-3D 
dimensional  crossover.  We  have  recently  detected  antiferromagnetic  order  in  the  3D  Fermi- 
Hubbard  model  using  spin-dependent  Bragg  scattering  of  light.  This  was  a  major  goal  of  our 
project,  enabled  by  the  development  of  a  novel  scheme  to  evaporatively  cool  atoms  in  an  optical 
lattice.  We  have  also  focused  on  developing  alternative  cooling  schemes.  We  have  also 
investigated  the  connection  between  quantum  integrability  and  thermalization  in  a  ID  Bose  gas, 
and  considered  equilibration  and  thermalization  phenomena  in  optical  lattices.  Other  theoretical 
issues,  such  as  the  effects  of  spin-orbit  coupling,  the  disordered  Bose-Hubbard  model,  and  many- 
body  localization,  have  also  been  explored. 


3D  Fermi-Hubbard  Model  at  Half-Filling 

The  Fermi-Hubbard  model  at  half-filling,  corresponding  to  one  particle  per  site,  is 
known  to  exhibit  antiferromagnetism.  Antiferromagnetism  is  a  property  shared  by  most  of  the 
undoped  parent  compounds  of  high-temperature  superconductors,  and  may  be  an  essential 
attribute  of  these  materials.  The  realization  of  antiferromagnetism  in  a  3D  optical  lattice, 
however,  has  been  impeded  by  the  extremely  low  temperatures  that  are  required.  The  Hulet 
group  has  developed  a  method  to  perform  evaporative  cooling  in  the  lattice  using  blue-detuned 
light  to  offset  the  confinement  of  the  red-detuned  optical  lattice.  In  this  scheme,  blue-detuned 
laser  beams  overlap  the  infrared  lattice  beams  along  the  three  orthogonal  directions,  as  shown  in 
Fig.  1.  Evaporation  can  be  controlled  by  tuning  the  chemical  potential  of  the  confined  atoms  to 


near  the  escape  threshold  with  the  intensity  of  the  blue-detuned  beams.  Significant  cooling  is 
observed.  In  addition,  the  blue-detuned  light  is  used  to  keep  the  density  fixed  at  half-filling 
during  the  lattice  ramp-up.  We  apply  this  scheme  to  a  3D  cubic  optical  lattice  filled  with  an 
equal  number  of  two  hyperfine  sublevels  of  the  ^Li  atom. 


Figure  1.  A  green  laser  beam  is  used  to  compensate  the  overall  confining  potential  from  a  conventional 
lattice  potential. 

The  beams  are  aligned  with  a  relative  precision  of  1  pm  using  the  three  identical  copies 
of  the  assembly  portrayed  in  Fig.  2.  The  purpose  of  the  compensating  beams  is  two-fold.  First, 
by  raising  the  chemical  potential  of  the  atoms  confined  to  the  lattice,  it  offers  the  opportunity  for 
atoms  to  evaporatively  cool  from  the  outer  regions  of  the  lattice,  as  depicted  in  Fig.  3.  Secondly, 
our  team  has  shown  theoretically  (Mathy,  Huse,  and  Hulet,  Phys.  Rev.  A  86,  023606,  2012)  that 
the  region  of  antiferromagnetism  may  be  enlarged  by  proper  choice  of  the  relative  intensities  and 
beam  sizes  of  the  green  beams,  relative  to  the  IR  beams.  Those  conditions  have  been 
incorporated  into  the  experiment. 
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In  a  red-detuned  3D  lattice:  In  a  compensated  lattice: 

(looking  along  a  lattice  direction) 


No  cooling 


Figure  3.  A  blue-detuned  compensating  beam  helps  to  improve  the  cooling  efficiency  inside  lattice. 

We  have  loaded  over  10^  atoms  into  the  lattice  at  a  magnetic  field  corresponding  to  a  scattering 
length  of  approximately  300  ao.  The  lattice  order  was  detected  by  Bragg  scattering  of  light  tuned 
to  near  the  2s-2p  resonance.  By  tuning  the  incident  angle  of  the  Bragg  beam  to  5 1°  with  respect 
to  a  lattice  plane  (001  direction)  and  collecting  the  scattered  light  with  a  lens,  we  were  able  to 
detect  a  huge  Bragg  scattering  signal.  The  signal  vs.  lattice  depth  fits  to  the  shape  given  by  the 
Debye-Waller  factor,  as  shown  in  Fig.  4.  This  experiment  demonstrates  that  Bragg  scattering 
produces  a  high-signal-to-noise,  easily  detectable  signal,  which  we  will  exploit  to  detect  anti¬ 
ferromagnetic  order. 

Bragg  scattering  from  the  symmetry  planes  corresponding  to  antiferromagnetic  order  is 
not  as  strong  or  as  forgiving  as  the  001  Bragg  scattering.  In  order  to  be  spin  sensitive  we  tune 
the  Bragg  probe  in  between  the  two  spin  states  corresponding  to  the  pseudo-spin  V2  system.  No 
signal  is  present  without  some  magnetic  ordering,  and  this  only  begins  when  the  system  is 
sufficiently  close  to  the  Neel  transition.  This  transition  temperature  is  predicted  by  theory  to  be  a 
maximum  of  -0.03  Tp.  By  adjusting  the  intensity  of  the  compensating  green  beams  to  be  as 
large  as  possible  without  significantly  distorting  the  overall  confinement  shape,  we  detect  a 
Bragg  signal  that  is  -50%  above  the  background  of  diffuse  scattering.  A  “rocking  curve” 
showing  the  angular  distribution  of  this  Bragg  signal  is  shown  in  Fig.  5.  The  angular  width 
corresponds  to  a  correlation  length  of  -4  sites,  which  quantum  Monte  Carlo  calculations  show  is 
a  temperature  that  is  perhaps  twice  T^.  By  adjusting  the  parameters  of  the  green  compensating 
beams  we  hope  to  cool  further  and  to  sharpen  the  Bragg  signal. 
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Figure  4.  Bragg  signal  from  the  3D  square  lattice  potential. 


Angle  along  chamber 


Figure  5.  Rocking  curve  (angular  distribution)  of  the  antiferromagnetic  Bragg  signal. 
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Trivedi’s  group  has  analyzed  the  magnetic  correlations  in  the  3D  Hubbard  Hamiltonian 
using  determinantal  quantum  Monte  Carlo  on  up  to  10  lattices,  including  the  effects  of 
population  imbalance.  Their  aim  is  to  guide  the  Rice  experimental  search  for 
antiferromagnetism  by  calculating  the  temperature  dependence  of  the  antiferromagnetic  structure 
factor  and  thereby  determining  the  temperature  scale  at  which  recent  experiments  might  be 
operating.  They  show  that  changing  the  density  away  from  half-filling,  or  the  spin  polarization 
away  from  zero,  within  the  expected  experimental  uncertainties,  does  not  significantly  affect  the 
results  at  the  identified  temperature.  By  comparing  simulations  of  the  Hubbard  and  Heisenberg 
models  they  have  determined  the  extent  of  finite  size  effects  close  to  the  Neel  transition.  Figures. 
6  and  7  plot  the  antiferromagnetic  structure  factor  as  functions  of  interaction  strength  U  and 
temperature  T,  respectively. 


Figure  6.  For  a  fixed  temperature  T,  the  antiferromagnetic  structure  factor  is  peaked  at 
U~8=10t  reflecting  the  non-monotonic  behavior  of  TNeei  which  peaks  at  this  optimal  U  value. 
This  suggests  the  optimal  value  that  should  be  used  in  experiments. 
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Figure  7.  Scaling  behavior  of  the  AF  structure  factor  when  plotted  as  a  function  of  T/TNeei  for 
various  U. 


Liu’s  group  has  investigated  physics  beyond  the  single -band  model.  They  discovered 
that  the  mixing  between  even  and  odd  orbitals,  such  as  s  and  p,  is  equivalent  to  a  partial  Rashba 
spin-orbital  coupling  when  the  orbital  space  (s,  p)  is  identified  as  pseudo-spin  1/2,  hence 
emergent  effective  spin-orbital  coupling.  They  further  extended  this  study  to  the  case  of  a  two- 
dimensional  array  of  such  ladders  (i.e.,  a  2D  array  of  ID  atomic  “tubes”).  Upon  crossing  over  to 
such  a  quasi-two-dimensional  system,  the  edge  modes  from  individual  ladder  form  a  parity- 
protected  flat  band  at  zero  energy.  Those  zero  modes  are  shown  to  be  the  equivalent  of 
decoupled  Majorana  fermions. 


One-Dimensional  Polarized  Fermi  Gas 

We  continue  to  upgrade  the  apparatus  used  to  produce  polarized  Fermi  gases,  and  the 
laser  system,  in  particular.  Our  focus  is  the  direct  detection  of  the  FFLO  superfluid  state.  Our 
approach  is  to  detect  the  non-zero  center-of-mass  momentum  pairs  that  are  the  hallmark  of  the 
FFLO  mechanism.  To  do  this,  we  have  constructed  a  blue-detuned  “anti-confining”  laser  system 
that  flattens  the  axial  harmonic  potential,  while  maintaining  confinement  radially.  This  laser 
system  is  nearing  completion,  and  we  expect  to  measure  time-of-fiight  momentum  distributions 
in  the  next  few  months. 


6 


Theoretical  studies  of  expansion  dynamics  of  ID  Fermi  gas  have  been  carried  out  by 
Bolech's  and  Pu's  groups  using  unbiased  numerical  techniques  such  as  density-matrix- 
renormalization-group  and  time-evolving-block-decimation  methods.  Such  studies  show  that  the 
time-of-flight  expansion  in  ID  is  highly  nontrivial  and  provide  better  understandings  of  the 
signature  of  the  FFLO  pairing. 

The  dimensional  crossover  is  another  important  phenomenon  that  we  are  investigating 
using  this  system.  The  phase  diagram  of  a  trapped  Fermi  gas  with  population  imbalance  is 
sensitive  to  the  dimensionality.  Liu's  group  tackled  this  problem  by  investigating  the  ID 
dynamical  evolution  of  an  initial  state  prepared  to  mimic  the  ground  state  of  the  3D  system. 
Bolech's  group  considered  the  coupling  between  an  array  of  ID  tubes  and  performed  a  mean- 
field  calculation  based  on  the  Bogoliubov-de  Gennes  equations  to  find  the  phase  diagram  of  the 
system.  It  is  shown  that  the  magnetization  of  the  system  undergoes  an  incompressible- 
compressible  transition  as  a  function  of  magnetic  field  and  intertube  tunneling  strength  (see  Fig. 
8). 


Figure  8.  (a)  Illustration  of  the  coupling  between  a  2D  array  of  tubes,  (b)  Filling  imbalance  per  tube,  M, 
in  a  radially  uniform  system  as  a  function  of  chemical  potential  difference  (or  magnetic  field)  ,  h,  and 
intertube  tunneling,  t.  The  lobe -like  plateaus  represent  magnetic  incompressible  states  with 
corresponding  integer  filling,  while  the  outside  region  has  finite  compressibility.  [From  Sun  and  Bolech, 
Phys.  Rev.  A  85,  051607(R)  (2012)] 
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Equilibration  and  Thermalization 

The  development  of  new  lattice  cooling  techniques  will  play  a  critical  role  in  experiments 
that  aim  to  probe  quantum  phases  that  are  currently  out  of  reach.  DeMarco’s  group  has 
developed  a  new  Raman  cooling  technique  that  works  for  any  atomic  species  trapped  in  an 
optical  lattice  and  does  not  require  a  spin-dependent  lattice  or  bath.  The  basic  scheme,  which  is 
closely  related  to  traditional  evaporative  cooling,  is  shown  in  Fig.  3a.  The  atoms  are  trapped  in  a 
lattice  and  are  prepared  in  one  hyperfine  ground  state  F.  Two-photon  stimulated  Raman 
transitions  are  used  to  transfer  high  quasi-momentum  atoms  to  an  excited  band  in  another 
hyperfine  state  F’.  Those  atoms  are  then  removed  using  near-resonant  light  that  only  affect 
atoms  in  state  F”.  The  remaining  atoms  collide  and  rethermalize  to  lower  temperature.  The 
cooling  rate  for  this  technique  can  be  very  high,  since  the  Raman  rate  is  limited  by  the  excited- 
state  bandwidth,  which  is  one  to  two  orders  of  magnitude  larger  than  that  of  the  ground  band. 


\F,mp) 

Figure  9.  Raman  cooling  method  shown  in  terms  of  the  band  structure.  The  general  technique 
(a)  involves  Raman  lasers  (light  red  arrows)  driving  atoms  between  hyperfine  states  in  a  quasi¬ 
momentum  selective  way.  Proof  of  principle  data  was  acquired  showing  the  variant  in  (b),  in 
which  microwaves  (gray  arrow)  are  used  to  transfer  atoms  between  hyperfine  states. 

DeMarco’s  group  demonstrated  this  technique  in  a  proof-of-principle  measurement  using 

07 

Rb  atoms  trapped  in  a  spin-dependent  lattice.  Figure  9b  shows  this  variant  of  the  method, 
which  involves  transferring  atoms  from  ground  to  the  first  excited  band  in  the  1 1,  -l)  state  using 
a  stimulated  Raman  transition.  Adiabatic  rapid  passage  driven  by  a  microwave-frequency 
magnetic  field  is  used  to  transfer  only  the  atoms  in  the  excited  band  to  the  1 2, 0)  state  (which 
does  not  experience  the  lattice),  where  they  are  removed  using  resonant  light. 
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Typical  images  before  and  after  a  single  cooling  step  are  shown  in  Fig.  10.  The  quasi¬ 
momentum  spread  of  the  atoms  is  reduced  and  condensed  fraction  is  increased.  Data 
demonstrating  this  cooling  effect  are  shown  in  Fig.  1 1  for  different  hold  times  in  the  lattice  after 
the  cooling  is  performed.  This  measurement  represents  the  first  demonstration  of  direct  cooling 
quasi-momentum  in  an  optical  lattice. 


Figure  10.  Images  showing  evidence  for  cooling.  Data  is  shown  for  a  gas  before  (a)  and  after 
(b)  a  cooling  pulse  is  applied.  In  (c),  traces  through  the  image  along  a  lattice  direction  (dashed 
line)  are  shown  before  (red)  and  after  (blue)  cooling. 


Figure  11.  Data  showing  cooling  for  different  hold  times  after  the  cooling  pulse.  Data  are 
shown  with  (blue)  and  without  (red)  a  cooling  pulse.  The  lines  indicate  the  average  and  standard 
deviation  taken  across  hold  times.  Cooling  increases  condensate  fraction  Nq/  N  and  reduces  the 

rms  size  of  the  gas,  which  is  proportional  to  yjr  . 

The  primary  limitation  on  this  method  is  sensitivity  to  magnetic  fields  at  the  mG  level. 
DeMarco’s  group  is  now  developing  a  laser  system  to  implement  the  most  general  version  of  this 
scheme  (Fig.  3a)  using  the  Raman-driven  magnetic-field- insensitive  transition  between  the  |l,0) 

and  1 2,0)  states.  A  sweep  of  the  Raman  detuning  will  be  used  to  continuously  cool  the  gas. 
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In  a  strong  theory-experiment  collaboration,  the  DeMarco  group  and  the  Mueller  group 
studied  the  dynamics  of  the  number  of  condensed  atoms,  following  a  lattice  ramp.  In  a  work 
published  in  Physical  Review  A  as  a  rapid  communication  [Stefan  S.  Natu,  David  C.  McKay, 
Brian  DeMarco,  and  Erich  J.  Mueller,  Phys.  Rev.  A  85,  061601],  they  showed  that  the 
condensate  fraction  evolves  significantly  on  a  100  microsecond  timescale.  A  positive 
consequence  of  this  rapid  response  is  that  local  equilibrium  is  readily  maintained  during  lattice 
ramps.  A  negative  consequence  is  that  "bandmapping"  techniques,  which  rely  on  changing  the 
lattice  parameters  quickly  compared  to  the  response  time,  are  not  able  to  give  quantitative 
information  about  such  systems.  These  results  will  be  fed  back  into  our  attempts  to  produce 
better  cooling  protocols.  Further  theoretical  support  came  from  Ceperley's  group  who,  using 
quantum  Monte  Carlo  techniques,  have  computed  entropies,  condensate  fractions  and 
momentum  distributions  for  bosons  in  3D  optical  lattices  and  are  completing  a  comparison  with 
DeMarco's  experiments.  This  will  be  used  to  better  characterize  experimental  densities  and 
temperatures. 

In  related  theoretical  studies  [Stefan  S.  Natu,  Erich  J.  Mueller,  arXiv:  1201.6674],  the 
Mueller  group  used  perturbation  theory  to  study  both  density-density  correlations  and 
momentum  distributions  of  a  Bose  gas  after  suddenly  turning  off  interactions  in  the  lattice.  They 
found  that  correlations  spread  ballistically,  with  a  speed  given  by  the  sharpest  part  of  the  single¬ 
particle  dispersion.  Quite  surprisingly,  they  found  that  the  ramps  to  zero  interactions  shared  all 
of  the  qualitative  features  of  experiments  in  the  strongly  interacting  regime  [Cheneau  et  al. 

Nature  481,  484  (2012)]  —  interactions  only  renormalized  the  propagation  speed.  They  are  now 
using  Bogoliubov  theory  to  push  these  calculations  to  stronger  interactions. 

In  collaborating  with  DeMarco,  Ceperley’s  group  have  used  quantum  MC  to  further 
studies  of  several  different  lattice  systems.  The  first  involves  comparisons  with  experiments 
conducted  at  DeMarco  group  for  the  clean  Bose-Hubbard  model  in  a  trap.  This  study  involved  N 
~  200,000  particles  and  focused  on  possible  thermometry  in  experiments  by  calibrating  against 
QMC  results.  The  work  has  been  completed  and  will  be  sent  for  publication  shortly;  it  calls  into 
question  the  equilibration  of  these  systems  at  high  entropy  near  .  Their  second  project  is  on 
disordered  systems.  They  have  performed  ab  initio  calculations  to  compare  directly  with 
experiment  -  the  first  of  its  kind  and  can  explain  the  absence  of  the  re-entrant  superfluid 
phenomena.  They  are  currently  in  the  process  of  understanding  the  different  features  arising  in 
the  trap  induced  heterogeneity,  e.g.,  how  the  different  domains  coexist  and  transform.  Towards 
this  end  they  have  undertaken  massive  parallel  computations  of  the  phase  diagram  including 
disorder  averaging. 

David  Huse's  group  is  interested  in  spin  and  energy  transport  in  universal  Fermi  gases. 
They  investigated  the  superdiffusive  nonequilibrium  dynamics  of  a  single  impurity  atom  in  a  low 
temperature  single-species  Fermi  gas  [Hyungwon  Kim,  David  A.  Huse,  Phys.  Rev.  A  85,  043603 
(2012)],  and  also  worked  out  the  diffusive  transport  behavior  of  the  polarized  Fermi  gas, 
including  heat  transport,  spin  Seebeck  and  spin  Peltier  effects.  In  addition,  they  also  studied 
nonequilibrium  critical  dynamics  of  the  quantum  Ising  spin  chain  and  of  the  Boson  Mott 
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insulator  in  a  tilted  potential  [Michael  Kolodrubetz,  Bryan  K.  Clark,  David  A.  Huse,  Phys.  Rev. 
Lett.  109,  015701  (2012)]  ,  showing  universality  between  the  cold-atom  "model"  studied  in 
Greiner's  lab  and  the  Ising  spin  chain. 

Huse’s  group  continued  their  study  on  many-body  localization.  They  have  found  a  new 
statistical  mechanics  of  phases  and  phase  transitions  within  the  many-body  localized  regime. 
These  phases  and  phase  transitions  show  up  in  the  dynamics  and  in  the  exact  many-body 
eigenstates  of  an  isolated  quantum  many-body  system,  but  are  invisible  to  the  equilibrium 
thermodynamics.  The  new  ordered  phases  (including  topologically  ordered  phases)  are 
dynamically  protected  by  the  many-body  localization  against  their  equilibrium  instability  to 
thermal  fluctuations.  This  opens  up  new  possibilities  for  stable  quantum  memories. 

In  another  work,  they  explained  a  sense  in  which  a  system  whose  eigenstates  are  all 
many-body  localized  can  be  viewed  as  an  integrable  system  with  localized  conserved  operators. 
These  localized  conserved  operators  do  interact,  and  this  causes  dephasing,  and  entanglement 
that  spreads  logarithmically  in  time. 

They  also  explored  eigenstate  phase  transitions  that  occur  within  the  ferromagnetic  phase 
of  an  infinite -range  Ising  model  that  obeys  the  Eigenstate  Therm ah/ati on  Hypothesis:  the  exact 
eigenstates  are  "Schrodinger  cat  states"  and  have  phase  transitions  in  the  dynamics  of  how  they 
tunnel  or  thermally  excite  themselves  to  flip  between  the  "up"  and  "down"  magnetized  phases. 
They  are  also  exploring  nonequilibrium  steady  states  with  energy  currents  flowing  between  two 
reservoirs. 

During  the  past  few  months,  Weiss’  group  has  finally  disentangled  all  heating  and  loss 
processes  in  their  Quantum  Newton's  Cradle  experiments,  and  they  found  a  way  to  calibrate  the 
expected  rate  of  diffractive  three-body  collisions,  by  referencing  it  to  the  3-body  inelastic  loss 
rate.  They  see  evidence  for  the  onset  of  thermalization  due  to  weak  non-integrability  in  a  ID  gas, 
but  with  a  time  constant  that  is  ~10  times  slower  than  conventional  theory  predicts  (see  Fig.  12). 
This  suggests  that  the  conventional  statistical  mechanics  that  underpins  the  Boltzmann  equation 
cannot  be  applied  too  close  to  integrability,  and  strongly  hints  that  there  may  be  a  quantum  KAM 
regime,  a  regime  away  from  integrability  where  systems  still  do  not  conventionally  thermalize. 
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Figure  12.  Evolution  of  the  total  energy  in  a  dense  quantum  Newton’s  cradle  [ref;  T.  Kinoshita, 
T.  Wenger,  and  D.S.  Weiss,  “A  quantum  Newton’s  cradle”.  Nature  440,  900  (2006)].  The  dip  at 
early  times  results  from  evaporative  cooling  due  to  diffractive  3-body  collisions,  while  the 
energy  increase  comes  mostly  from  spontaneous  emission.  The  cooling  rate  is  10  times  slower 
than  predicted  for  these  conditions.  Because  of  the  absence  of  diffractive  2-body  collisions  in 
this  system,  this  may  be  also  the  first  observation  of  elastic  3 -body  collisions. 


Other  Theoretical  Issues 

In  addition  to  those  mentioned  above,  our  team  have  explored  various  theoretical  issues 
that  are  not  necessarily  tied  closely  to  the  experiments  carried  out  in  the  Hulet  and  DeMarco  labs, 
but  are  broadly  related  to  the  OLE  program. 

Ho's  group  has  continued  the  study  of  strongly  interacting  Fermi  gases,  with  focus  on 
obtaining  thermodynamic  information  of  strongly  interacting  Fermi  gas  through  its  density 
profde.  In  the  previous  funding  periods,  Ho's  group  has  discovered  algorithms  to  process  the 
density  to  obtain  thermodynamic  functions  and  to  determine  phase  transitions  at  finite 
temperatures  and  zero  temperature  (7=0).  His  group  has  demonstrated  the  algorithm  for 
attractive  ID  Fermi  gas,  which  has  only  7=0  but  not  finite  7  transitions.  In  the  last  funding 
cycle,  they  have  extended  this  study  to  repulsive  Fermi  gas,  which  is  related  to  the  controversy  of 
the  discovery  of  ferromagnetism  reported  by  the  Ketterle  group.  They  have  shown  that  all  the 
purported  phenomena  for  ferromagnetism  can  be  accounted  for  by  Pauli  blocking  physics  in 
the  so-called  upper  branch  of  the  energy  spectrum.  In  this  funding  cycle,  they  have  extended  this 
study  to  Fermi  gases  with  narrow  resonances  and  to  repulsive  Bose  gases.  In  addition,  they  have 
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shown  recently  that  ferromagnetism  in  a  spin- 1/2  Fermi  gas  is  the  (equal  mass)  limiting  case  of 
phase  separation  in  a  binary  mixture  of  Fermi  gases,  and  have  proven  a  theorem  that  phase 
separation  (the  analog  of  ferromagnetism)  must  occur  for  sufficiently  large  mass  asymmetry.  In 
addition,  Ho's  group  continues  to  explore  new  ways  to  increase  the  energy  scale  of  lattice  gas 
systems  by  exploring  the  interaction  between  a  lattice  gas  with  another  atomic  species  not 
confined  by  the  lattice  which  acts  as  an  agent  that  induces  effective  long-range  interaction. 

Ceperley's  group  has  done  exact  numerical  calculations  of  few  fermion  systems  where 
scattering  length  is  on  the  order  of  the  inter-particle  spacing,  which  is  important  since  many  of 
the  OLE  experiments  are  done  in  this  regime.  They  have  seen  how  underlying  molecular  states 
influence  scattering  states  and  determined  that  the  ferromagnetic  transition  density  is  changed  by 
a  factor  of  two  if  one  includes  the  underlying  molecular  states,  in  agreement  with  the  experiment 
of  Jo  et  al.  They  have  developed  techniques  to  generalize  this  technique  from  the  lattice  model 
into  the  continuum. 

Liu's  group  has  investigated  topological  phases  of  fermions  in  optical  lattices.  This 
research  is  broadly  related  to  the  Fermi-Hubbard  model  project,  which  is  part  of  the  core  thrusts 
of  the  OLE  program.  Cold  fermionic  atoms  were  discovered  to  exhibit  topological  semimetal 
and  insulator  phases,  both  being  related  to  the  current  on-going  hot  topics  in  Condensed  Matter 
physics. 

Both  Mueller's  and  Pu's  groups  have  explored  the  effects  of  non- Abelian  synthetic  Gauge 
field  on  cold  atoms.  Mueller's  group  has  conducted  a  number  of  numerical  experiments  to 
evaluate  the  feasibility  of  observing  non-Abelian  statistics  with  cold  atoms.  In  a  paper  published 
in  Physical  Review  Letters  [Eliot  Kapit,  Paul  Ginsparg,  and  Erich  Mueller,  Phys.  Rev.  Lett.  108, 
066802  (2012)]  they  numerically  implemented  a  long-discussed  thought  experiment.  In  the 
experiment  one  imagines  having  a  system  in  the  fractional  quantum  Hall  regime  (an  exotic  state 
found  in  semiconductor  heterostructures  in  large  magnetic  fields).  One  uses  two  lasers  to 
produce  two  potential  wells.  One  then  imagines  slowly  moving  the  laser  spots,  switching  the 
location  of  the  wells.  According  to  variational  calculations,  one  should  be  able  to  perform 
quantum  computation  by  a  sequence  of  such  moves.  Kapit  et  al.  numerically  conducted  this 
experiment,  verifying  that  it  works.  They  also  proposed  a  cold-atom  experiment  which  would 
demonstrate  this  physics.  Pu's  group  has  studied  the  effects  of  spin-orbit  coupling  in  quantum 
gases  induced  by  the  non-Abelian  Gauge  field.  They  have  shown  that  a  properly  engineered 
spin-orbit  coupling  term  can  greatly  enhance  fermionic  pairing  and  may  lead  to  exotic 
topological  phases  in  both  Boson  and  Eermi  gases. 

Both  Bolech  and  Trivedi  studied  Bose-Hubbard  model.  Bolech  studied  a  Bose-Hubbard 
model  having  on-site  repulsion,  nearest  neighbor  tunneling,  and  ferromagnetic -like  coupling 
between  the  nearest-neighbors'  occupation  parities.  Eor  a  uniform  system  in  any  dimension  at 
zero  tunneling,  they  obtained  an  exact  phase  diagram  characterized  by  Mott  insulator  (MI)  and 
pair  liquid  phases  and  regions  of  phase  separation  of  two  Mis.  Eor  a  general  trapped  system  in 
one  and  two  dimensions,  they  performed  quantum  Monte  Carlo  and  Gutzwiller  mean-field 
calculations.  Both  of  which  show  the  evolution  of  the  system  from  a  superfluid  to  a  wedding- 
cake  structure  of  Mis  with  the  occupations  jumping  by  2,  as  the  parity  coupling  increases.  They 
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also  identified  an  exotic  pair  superfluid  at  a  relatively  large  tunneling  strength.  Their  model 
ought  to  effectively  describe  recent  findings  in  imbalanced  Fermi  gases  in  two  dimensional 
optical  lattices. 

Motivated  by  the  experimental  realization  of  synthetic  spin-orbit  coupling  for  ultracold 
atoms,  Trivedi’s  group  have  investigated  the  phase  diagram  of  the  Bose-Hubbard  model  in  a 
non- Abelian  gauge  field  in  two  dimensions  (a  is  the  spin-orbit  coupling  strength  and  A,  is  the 
relative  strength  of  bosons  of  opposite  spin).  Using  a  strong  coupling  expansion  in  the  combined 
presence  of  spin-orbit  coupling  and  tunable  interactions,  they  find  a  variety  of  interesting 
magnetic  Hamiltonians  in  the  Mott  insulator  (MI),  which  support  magnetic  textures  such  as  spin 
spirals  and  vortex  and  Skyrmion  crystals.  An  inhomogeneous  mean-field  treatment  shows  that 
the  superfluid  (SF)  phases  inherit  these  exotic  magnetic  orders  from  the  MI  and  display,  in 
addition,  unusual  modulated  current  patterns.  They  present  a  slave-boson  theory  which  gives 
insight  into  such  intertwined  spin-charge  orders  in  the  SF,  and  discuss  signatures  of  these  orders 
in  Bragg  scattering,  in  situ  microscopy,  and  dynamic  quench  experiments. 

Several  groups  have  investigated  the  properties  of  quantum  gases  subject  to  synthetic 
spin-orbit  coupling.  Together  with  his  collaborators  (Ben  Lev,  Xiaoling  Cui,  and  Hui  Zhai),  Ho 
has  shown  that  the  heating  due  to  the  added  laser  fields  which  generate  the  spin-orbit  coupling  is 
substantially  reduced  if  one  uses  Lanthanide  atoms  rather  than  the  alkalis.  Both  Mueller’s  and 
Pu’s  group  have  studied  the  “Majorana  modes”  in  spin-orbit  coupled  topological  Fermi 
superfluids.  Ho’s  group  also  completed  a  work  on  the  spin-orbit  effects  in  ID  repulsive  Fermi 
gas.  It  was  shown  that  at  infinite  repulsion,  where  the  system  has  infinite  spin  degeneracy,  even 
a  weak  spin-orbit  interaction  can  lead  to  the  generation  of  exotic  phases  such  as  a  large  spin 
spiral. 
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